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ABSTRACT 

We present panoramic it' and optical ground-based imaging observations of a complete 
sample of low- redshift (0 < z < 0.4) galaxies selected in the ultraviolet (UV) at 
> 1 2000 A using the balloon-borne FOCA instrument of Milliard et al. This survey is 

highly sensitive to newly-formed massive stars, and hence to actively star-forming 
galaxies. We use the new data to further investigate the optical, stellar population 
and star formation properties of this unique sample, deriving accurate galaxy types 
and /c-corrections based on the broad-band spectral energy distributions. 

When combined with our earlier spectroscopic surveys, these new data allow us 
\ to compare star-formation measures derived from aperture-corrected Ha line fluxes, 

UV(2000A) andu'(3600A) continuum fluxes on a galaxy- by-galaxy basis. As expected 
' from our earlier studies, we find broad correlations over several decades in luminos- 

ity between the different dust-corrected star-formation diagnostics, though the scatter 
Q ■ is larger than that from observational errors, with significant offsets from trends ex- 

| pected according to simple models of the star-formation histories (SFHs) of galaxies. 

c/3 . Popular galaxy spectral synthesis models with varying metallicities and/or initial mass 

functions seem unable to explain the observed discrepancies. 

We investigate the star-formation properties further by modelling the observed 
spectroscopic and photometric properties of the galaxies in our survey. We find that 
' nearly half of the galaxies surveyed possess features that appear incompatible with 

^ \ simple constant or smoothly declining SFHs, favouring instead irregular or temporally- 

varying SFHs. We demonstrate how this can reconcile the majority of our observations, 
enabling us to determine empirical corrections that can be used to calculate intrinsic 
star formation rates (as derived from Ha luminosities) from measures based on UV (or 
u') continuum observations alone. We discuss the broader implications of our finding 
that a significant fraction of star-forming galaxies have complex SFHs, particularly in 
the context of recent determinations of the cosmic SFH. 

Key words: surveys - galaxies: evolution - galaxies: luminosity function, mass func- 
tion - galaxies: starburst - cosmology: observations - ultraviolet: galaxies 



1 INTRODUCTION 

An accurate determination of the star-formation history 
(SFH) of the Universe is one of the key goals of modern ob- 
servational cosmology. Studies which measure, analyse and 
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model its precise form are important not only in indicating 
likely epochs of dominant activity, but also for comparisons 
with the predictions of semi-analytical models of galaxy for- 
mation, and as such have provided a major impetus towards 
a fuller understanding of the physical mechanisms of galaxy 
evolution. 

Various techniques now exist to measure star-formation 
rates (SFRs) in distant galaxies, all of which are sensi- 
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tive in some way to the number of young, short-lived, 
and hence massive stars. These include nebular recombi- 
nation (e.g. Ha) and forbidden line ( e.g. [Qui) emission 
jKennicuttlll983l: [Gallagher et al.lll989HGallego et alJll995t 
iTresse et alJ 120021: Iffippelein et alJl2003h . ultraviolet (UV) 



jDonas et alJ Il987l 


iLillv et alJ 


19961 IConnollv et al 


1997; 


Trever et alJ 119981 


Cowie et al. 


19991: ISullivan et al. 


2000; 


Wilson et al.l 120021 


, decimetric radio emission JCondod 


19921: ICram et alJ 


Il998l: lHaarsma et al.l l2000h. and far- 



IBlain et aflll999l:lKe"wlev et al.ll2002l) (for a revi ew of these 
techn iques and their associated uncertainties, see iKennicutd 
Currently, many of these diagnostics are only avail- 



1199; 

able over limited (and sometimes non-overlapping) red- 
shift ranges, so studies of the evolution of the cosmic SFH 
have tended to combine many disparate studies measuring 
star-formation using a selection of the techniques described 
above, without an accurate understanding of how underlying 
physical processes in galaxies can affect different diagnostics. 

However, there has been recent progress in inter- 
comparing and calibrating these measures by investi- 
gating the degree to which the different diagnostics 
agree when measured for the s ame galaxies. Work to 
date comparing Hq/UV measures jGlazebrook et alJll999t 



date comparing Ha/UV measures (Cclazcbrook et al. 1999; 
Sullivan et alJl2000t Irfopkins et aliboOlt iBell fc Kennicutd 



20011; iBuat et all 120021) and radio/Ha/UV measures 
JCram et alll998tlHopkins et aljl200 it ISullivan et alJl200lTl 

has revealed broad correlations but with a significant scat- 
ter, together with some offsets from relations expected from 
simple SFH scenarios. Such discrepancies indicate either 
poorly understood dust extinction corrections (particularly 
on UV measures), some physical parameter varying from 
galaxy to galaxy, such as the metallicity or initial mass func- 
tion, or variations in the timescale of recent star-formation 
and consequent contamination of UV-continuum derived 
SFRs by older and less massive stars. 

We have been using a unique sample of local UV- 
selected galaxies to investigate these issues, construc ted us- 
ing the balloon-borne camera of Milliard^tjalJ_fcl92) which 
images the sky at 2000 A. In iTrever et al] (Il99gfi we used 
this sample to construct the first local UV luminosity func- 
tion, finding an integrated star-formation density higher 
than t hat found in local emission-line surveys. ISullivan et al] 
( 2000) extended this sample, and performed initial investi- 
gations into the physical nature of star-formation in the sam- 
ple comparing Hq and UV-derived SFRs. Evidence for non- 
linearities and significant scatter was found, though ques- 
tions over aperture corrections on the Hq measures and 
fc-corrections on the UV luminos ities remained. The ra- 
dio properties ijSullivan et aljl20o"lT) and chemical properties 
IContini et al]l2002ft of this sample have also been investi- 
gated. 

In this paper, we further constrain the nature of star- 
formation in the sample. We extend our previous studies to 
include new broad-band u' and optical Br 1 imaging observa- 
tions of a substantial fraction of the galaxies in the survey, 
allowing accurate aperture and fc-corrections to be made. 
We then use the combined samples to investigate the opti- 
cal and stellar population properties of the sample, as well 
as to investigate SFRs derived from a second UV contin- 



uum measure in conjunction with that previously measured 
at 2000 A. 

An outline of the paper is as follows. In Section |2] we 
introduce the redshift survey data, including the new imag- 
ing campaigns. Section [3] discusses the optical properties of 
the sample, allowing improved aperture and fc-corrections 
for our spectroscopic data. Section^Jthen examines the star- 
formation properties of the galaxies derived from Hq/UV/u' 
observations. We model and discuss the implications of our 
results in Section |H| and conclude in Section |S| Through- 
out this paper we assume a Qa = 0.7, Qm = 0.3, h — 0.70 
(where Ho = 100 /ikms -1 Mpc ) cosmology. 



2 THE SURVEY DATASETS 

This paper extends our earlier work of [ Trever et al] (Il998l . 
hereafter T98) and ISullivan et all feood hereafter S2000) 
by augmenting our photometric dataset with with more ex- 
tensive and precise optical photometry. Here, we introduce 
the new datasets which are analysed in this paper. We com- 
mence with a brief discussion of the existing datasets - UV 
imaging data taken from the FOCA experiment and the 
follow-up optical spectroscopy - followed by the details of 
the new imaging surveys of these fields. 

2.1 The FOCA redshift survey 

The FOCA instrument is a balloon-borne 40-cm Cassegrain 
with a single filter a pproximating a Gau ssian centred at 
2015A, FWHM 188A tMilliard et alJll992ft . Two fields ob- 
served by FOCA are studied here: the high galactic-latitude 
region of Selected Area 57 (SA57; RA = 13 h 06 m 10!43, 
Dec. = +29°0l'46'.'l, J2000) and the field of the z = 
0.0215 Leo cluster, Abell 1367 (A1367; RA = ll h 45 m 20:92, 
Dec. = +19°53'23'.'0, J2000). The UV exposure times cor- 
responded to a limiting magnitude of muv = 18.5 in the 
FOCA photometric system (m^y ~ 20.75). 1 Both fields 
have extensive multi-fibre optical spectroscopy obtained us- 
ing the WIYN/Hydra (3500 -6600 A; 3.1-arcsec diameter 
fibres) and WHT/WYFFOS (3500 -9000 A; 2.7-arcsec di- 
ameter fibres) telescope/instrument combinations, provid- 
ing redshifts for 224 UV-selected emission-line galaxies, of 
which Hq fluxes can be reliably measured in 111 objects. 
The wavelength coverage of the WIYN data preclude obser- 
vation of Hq. Fluxes, equivalent widths (EWs), and errors 
for each of the principle emission lines, as well as the strengt h 
of the 4000 A Balmer break (D4000; see e.g. lBruzuaJll983ft . 
are measured where the signal-to-noise (S/N) of the spectra 
permits. Where Hq is not detected, an upper flux limit is 
estimated at the emission wavelength using the noise char- 
acteristics of each spectrum. 

The Balmer emission lines are corrected for the effects 
of stellar absorption as described in earlier papers (S2000), 
with a median correction of 2.5 A. Fluxes are also corrected 
for the effects of Galactic extinction using the dust maps of 

1 Subsequent to the publication of S2000, the FOCA instrument 
team have re-calibrated the conversion required to transform 
FOCA instrumental magnitudes to an external absolute magni- 
tude system, resulting in a minor adjustment in the sense that 

TOncw = "1 S 2000 + 0- 02 - 
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ISchleeel et"afl (Il998l) and a lCardelli et alJ il989T> extinction 
law. These corrections are small, with typical E(B — V)s of 
0.01 and 0.02 for SA57 and A1367 respectively. The ratio of 
these corrected Ha and H/3 fluxes is used to determine the 
colour excess of the ionised gas, E(B — V) s , and hence the 
i nternal dust extinctio n assuming case-B recombination and 
a lCardelli et alJ ill 9891) attenuation law. AGN and QSO-like 
objects are removed from the star-forming sample on the ba- 
sis of their optical spectra (objects with broad emission lines 
are discarded). Of the remaining narrow-band objects, the 
vast majority are star-forming galaxies, as evidenced from 
diagnostic diagrams based upon the r atios of their princi- 
ple emission lines llContini et alj|2 002h. De tails of all these 
procedures can be found in lsTuhvanl i2002T) . 



2.2 Optical Imaging Campaigns 

The new optical data are taken from three sources: A w'/r' 
Palomar/LFC imaging programme of SA57, B-band imag- 
ing using the CFHT/12k camera, and supplementary gri 
data from the Palomar Digital Sky Survey (DPOSS) for both 
SA57 and A1367. 



2.2.1 The Palomar Large Format Camera 

The Large Format Camera ("LFC lSimcoe et al. 2000) on the 
Palomar 200-in Hale telescope is a mosaic of six 2048x4096 
pixel CCDs (though only four were available when the data 
for this paper were collected), covering a region of diame- 
ter ~ 24' with a pixel scale of 0.175 arcseconds. Data were 
collected over the course of six dark nights split over two 
observing runs in April 2000 and March 2001. The 1.5° di- 
ameter field of SA57 was mosaiced using Sloan Digital Sky 
Survey (SDSS) u' and r' filters, with a tessellation pattern 
avoiding bright foreground stars which created serious sat- 
uration and bleed trails on the LFC CCDs. Exposure times 
were 2400s in u' and 600s in r'. Due to time constraints, 
we concentrated our observations on the centre of the SA57 
field. 

The LFC data were reduced using a hybrid of the Na- 
tional Optical Astronomy Observat ories (NOAO ) mosaic re- 
duction software for use in iraf as well as 
custom-written routines to deal with some issues specific 
to the LFC. The data were bias-subtracted using the over- 
scan regions on the chips and nightly master bias frames. 
Flat-fielding was performed using dome flats (for r' data) 
and twilight sky flats (for u' data) taken on each night. We 
also constructed a 'super-sky' flat for the r' data by stack- 
ing all the r' observations; however the resulting sky flat- 
field frame was virtually indistinguishable from the dome 
flat, and hence was not used. The iraf routine mscskysub 
was used to remove any second-order sky gradients across 
the field by computing medians in boxes of size 50 pixels 
and fitting a low-order two dimensional function to the me- 
dian points. The residual 'fitted-sky minus sky mean' is sub- 
tracted from each pixel. 

We derived astrometric solutions for each chip in the 
two filters using astrometric calibration exposures. SA57 is 
a (deliberately chosen) high galactic latitude field where the 
surface density of bright foreground galactic stars required 
for astrometric calibration is low. We therefore observed as- 



trometric calibration fields taken from t he Astromet r ic Cal - 
ibration Regions (ACRs) catalogue of IStone et alJ il999t) . 
with star posi tions typically accurate to ±26 mas. We use 
Sextractor feertin fc Arnoutslll996D and the W CSTools 
software suite written by D. Mink (e.g. lMinkiri999l) to cre- 
ate approximate world co-ordinate system (WCS) solutions 
for each astrometric region in each filter based on stars in 
the ACR catalogues, discarding stars with large proper mo- 
tions. The WCS was then refined using the iraf program 
CCMAP and a TNX projection with a high order polyno- 
mial and half cross-terms. This flexible function was needed 
to fit the edges of the LFC chips where distortion can be 
significant. The plate solution for each filter/chip was then 
applied to every science exposure using the iraf mscred 
tasks, shifting the reference point for the plate solution for 
each science field using the USNO-A2.0 stars (i.e. there are 
enough of these stars to calculate a reference point shift, 
but not to derive a full plate solution). The final solution 
for each science field is accurate to 0.15-0.2", even for the 
regions of the chips that suffer significant distortion. 

Each exposure was projected and re-sampled onto a lin- 
ear WCS, and the individual dither steps at each pointing 
median-combined after masking cosmic-rays and chip de- 
fects to create one final dithered image for each chip, ap- 
plying multiplicative scaling factors to allow for airmass dif- 
ferences between dithers. The result is an astrometrically 
calibrated and photometrically stable combined science im- 
age. 

At the time the observing program was carried out, 
the SDSS filters were 'non-standard', in the sense that few 
published calibration stars were available. To avoid com- 
plex colour terms needed to place our data on the stan- 
dard a-Lyrae photometric system, we used standard spec- 
trophotometric stars to calibrate our data. On each night, 
two spectrophotometric standard stars were monitored at 
approximately two hour periods throughout the night to 
derive airmass-dependent extinction correction terms, aug- 
mented by dawn and dusk twilight observations of 4 - 5 other 
spectrophotometric standard stars. Our preferred standards 
were Feige 67 and Feige 34. We calculated the AB magni- 
tude of each standard star in the u' and r' filters using the 
spectrophot ometric star sp e ctra e nergy distributions (S EDs) 
taken from iMassev et alJ il988l) and iMassev fc Gronwalll 
Jl990f) . We convolved these SEDs with the filter and LFC- 
CCD responses (C. Steidel, private communication) to cal- 
culate the magnitude in a particular filter. 

The r' zero-point was extremely stable with a mean 
r.m.s. dispersion of ~ 0.01 mag. The u' zero-point was less 
well-defined, with an r.m.s. variation from standard star to 
standard star of ~ 0.03 mag. Due to the accuracy limit to 
which the spectra of these standard stars are calibrated an 
accuracy grea ter than 0.03 mag in u' is p robably not achiev- 
able (see e.g. ISteidel fc Hamilton! I1993T) . As an additional 
consistency check, the overlapping nature of the mosaicing 
technique allows us to compare magnitudes of objects in dif- 
ferent fields, and a night-to-night (or even year-to-year) com- 
parison is possible. We typically find agreement to within 
0.02 mag in r\ rising to 0.04 mag for the u' data. No offsets 
were seen between the same objects measured on different 
nights. 

Object catalog ues were created using SEXTRACTOR ver- 
sion 2.2.2 (IBertin fc Arnoutall996l) . We used the automatic 
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aperture magnitudes determined by SEXTRACT OR, similar 
to Kron's 'first moment' algorithm (lKronlll98Ctl . Our high 
galactic-latitude fields contain a low source-density of ob- 
jects (particularly in it'), so crowding on the fields was not 
a problem and de-blending was rarely required. 



2.2.2 CFHT and DPOSS data 

A second optical survey of SA57 was conducted using the 
Canada France Hawaii Telescope (CFHT) and the CFH12k 
instrument, a la rge mosaicing camera w ith an imaging area 
of 42'x28' fsee ICuillandre et ahlfeoOOl) . The observations 
were carried out on the 25th (3/4 night) and 26th (1/4 night) 
May 2000 using seven CFHT12k pointings to cover the SA57 
field. The first night was photometric, but poor weather on 
the second night curtailed the prog ram. Each poin ting com- 
prised two 300s dithered exposures. iLandoitl jl992T) standard 
star fields were interspersed with the science observations to 
achieve a photometric calibration of the data. 

The data were reduced at the TERAPIX centre using 
the TERAPIX 2 pipeline reduction software designed for use 
with the CFH12k camera according to a procedure described 
in the June 2002 TERAPIX progress report 3 . An image cat- 
alogue covering SA57 was constructed using the SEXTRAC- 
TOR package and here again we use the automatic aperture 
magnitudes. The catalogue is complete to B = 24. 

Finally, we augment the two dedicated observing pro- 
grams of SA57 with data taken from the Digitised Palomar 
Sky Survey (DPOSS). This provides d ata for both SA5 7 and 
A1367 in g, r and i Gunn filters fsee lGal et afll????ft . The 
areal coverage in both fields is 100 per cent, though the data 
are shallower (g < 22.5) and of limited photometric accuracy 
when compared to the Palomar/CFHT data, hence we use 
the DPOSS data where no other imaging is available. The 
optical sources in the various bands are uniquely matched 
within a search radius of 1.5". As in our earlier papers, we 
use a search radius of 10" to assign optical counterparts 
from our combined optical data to the FOCA sources. The 
problem of multiple optical counterparts is addressed in Sec- 
tion |XT] 



2.3 Comparing the magnitude systems 

To consistently inter-compare these new photometric 
data, we place all the magnitudes onto a system with 
the same calibration zeropoint. Our imaging data com- 
prises four different calibration systems: The FOCA sys- 
tem (e .g. |Milhaj^j^t_alJ ll992| ), the Vega system (e.g. 
IJohnson fc Morganlll953l : lFukugita et alJll995t). the AB sys- 
tem ( e.g. lOkell974) . and the Gunn system fe.g. lOke fc Gunn] 
ll98St) . We compute conversions between the different sys- 
tems using the SEDs of Vega, BD+17°4708 (for Gunn mag- 
nitudes), and appropriate filter and CCD response curves. 
Because of the advantages of a system where magnitudes are 
easy to interpret physically, we align all of our magnitudes 
onto the AB system, defined as: 



Table 1. Offsets used to convert the magnitudes to the AB sys- 
tem. 



Filter 


(A) 


Conversion from 
foca — > AB vega — * AB gunn — > AB 


FOCA-UV 


2015 


+2.25 


Palomar-u' 


3603 


+0.92 


CFHT-B 


4407 


-0.10 


DPOSS g 


5330 


+0.02 


Palomar-r' 


6282 


+0.17 


DPOSS r 


6908 


-0.19 



-2.51og(F„) -48.60 



2 http://terapix.iap.fr/ 

3 http:/ /terapix. iap.fr/doc/doc. html 



(1) 



where F is the flux in ergs -1 cm" 2 Hz -1 . For convenience, 
we list the conversion values between the various original 
systems to the AB system in Tabled 



3 OPTICAL PROPERTIES 

In this section we examine the photometric properties of the 
galaxy sample. We first describe our procedures for assign- 
ing optical counterparts to the FOCA source list, and then 
correct all of the observed magnitudes for the effect of inter- 
nal dust extinction. The new optical coverage also allows us 
to estimate aperture corrections for our spectral measures, 
as well as a cross-check on the accuracy of the flux calibra- 
tion of our spectra. We then continue with an exploration 
of the UV-optical properties, using the new data to derive 
improved galaxy types, fc-corrections and luminosities com- 
pared to our earlier analyses. 



3.1 FOCA — optical counterparts 

One question remaining from our earlier papers concerned 
the number of FOCA detections with no apparent counter- 
part in optical data, which would imply either blue UV — B 
colours or a number of false detections in the FOCA dataset. 
The optical imaging now available to us allows us to address 
this issue in a more comprehensive manner, as the depth of 
our B-band data is ~ 2 magnitudes deeper than the opti- 
cal dataset used in S2000. This magnitude limit enables us 
to detect all the optical counterparts to the FOCA source 
list. We cross-correlate the combined optical u Br' catalogue 
with the FOCA source catalogue using a search radius of 10" 
(see e.g. T98 or S2000 for details of the astrometric preci- 
sion of the FOCA data). For each FOCA object, we record 
the number of matches, and assign each the nearest optical 
detection. 

The number of multiple counterpart cases (where more 
than one optical detection lies within the search radius) is 
increased compared to S2000; however the number of UV 
detections with no optical counterpart is smaller: 50 per 
cent of the FOCA sources have multiple counterparts whilst 
less than 6 per cent have none. This 6 per cent thus rep- 
resents a reasonable lower limit to the false detection rate 
in the FOCA source list, as the magnitude limit of the new 
data allows us to detect FOCA counterparts ~ 1 magnitude 
bluer in UV — B than the bluest unreddened model starburst 
galaxy SED. The extra number of counterparts we find here 
are therefore most likely chance associations resulting from 
the higher surface density of objects in the fields: the bright- 
est optical counterpart is also the closest in 85 per cent of 
the multiple counterpart cases. Where the brightest optical 
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detection is more than 3 magnitudes brighter then the next 
brightest counterpart, we assume that the contribution to 
the UV flux of the fainter object is negligible and ignore 
it in our analyses. In our following analyses, 'multiple op- 
tical counterparts' galaxies are defined as those FOCA ob- 
jects with two optical counterparts within B-band measures 
within 3 magnitudes of each other. 



3.2 Dust and fc-corrections 

There are a variety of extinction laws (e.g., Milky Way, 
LMC, SMC) that could be used to apply the dust correc- 
tions factors derived from the Balmer decrements to the UV 
and optical magnitudes, giving quite different extinction cor- 
rections, particularly at UV wavelengths. The reddening of 
the UV continuum depends sensitively on the geometrical 
details of the dust-star-gas mix, with the reddening of the 
stellar continuum likely to be different to the obscuration of 
the ionised gas (derived empirically from the nebular emis- 
sion via the Balmer decrement), as the stars and gas may 
occupy different areas wi thin a galaxy with differing dust 
cover i ng factors (see e.g . [Fuiel li et alJ fl988t ICalzetti et alJ 
11994 iMas-Hesse fc Kunthlll999t) . As in our previous work, 
we use the lclluzet^i^tHuTll200oh prescription, derived em- 
pirically from studies of local starburst gal axies, to 'cor- 
rect' our optical fluxes. ICalzetti et al.l i2000t) find that the 
colour excess of the stellar continuum, E(B — V) s , is re- 
lated to the colour excess of the ionised gas, E(B — V) e , as 
E(B - V) s = O.UE(B - V) s . 

We apply this reddening prescription to every filter 
passband in use in our survey using E(B — V) s derived from 
the Balmer decrement. Where the H/3 line is not detected, we 
use the SFR-depend ent correction derived for this sample of 
galaxies presented in lSullivan et al.l i200ll) using an empiri- 
cal SFR derived from the Ha line (where both the Ha and 
H/9 lines are absent, we use the median E(B — U) g =0.3 repre- 
sentative of the whole sample) . It is becoming apparent that 
luminosity-dependent dust corrections of this nature are not 
valid for star-formation across all Hubble types. However, for 
the mainly starburst or strongly star-forming galaxies typi- 
cal in this sample, the luminosity de pendent correction still 
appears an excellent approximation dAfonso et aljl2003T) . 

We convert our observed (and dust-corrected) magni- 
tudes, m, in the various passbands into rest-frame magni- 
tudes by writing the corrected absolute magnitude, M, at 
2 = as 

M = m obs - 5 log(^j^) - k(z) (2) 

where Dl is the luminosity distance at redshift z, m bs is the 
observed apparen t magnitude, and k(z) is the fc-correctio n 
(see, for example, lYoshii fc Takaharall98d lPoggiantilll997ft . 
We neglect any redshift-evolution in the physical properties 
of the galaxies in this calculation. 

The magnitude of the fc-corrections for template galaxy 
spectra using the range of the filters in this study strongly 
depends on galaxy type, with early-type galaxies (ellipti- 
cals and SOs) possessing larger fc-corrections in the optical 
passbands than galaxies with flatter spectra. However, fc- 
corrections in the FOCA passband are small with the largest 
fc-corrections arising from elliptical and SO type spectra, un- 



likely to be the predominant galaxy type in this UV-selected 
sample of objects. 

We assign each object in our survey a fc-correction 
by fitting the observed magnitudes to a set of model 
galaxy templates of varying type. Though we could do 
this by correlating our observed spectra with a set of tem- 
plate spectra, either via cross-correlation techniques (e.g . 
IConnollv et~al1ll995l IZaritskv et al]ll995t iHevl et al.lll997l) 
or using a principal compo nent analysis (see iFolkes et alJ 
ll999tlMadgwick et all2002l and references therein), this ap- 
proach requires spectra of a high S/N and wide wavelength 
coverage, which is not always the case in the spectra of our 
fainter objects. Instead, we use the popular technique of fit- 
ting a series of template SEDs to the broadband colours 
available for each galaxy and use the best-fitting SED for 
the estimation of fc-corrections in each filter. This was also 
the approach used in T98 and S2000 though based on only 
a single UV — b colour and using POSS-APM magnitudes, 
which have proved unreliable for bright galaxies. Here, we 
perform a more accurate calculation with the improved pho- 
tometric coverage. 

We determine our fc-corrections by fitting each colour 
for a given galaxy to SED set, linearly interpolating be- 
tween the SED types. These fits are combined by weight- 
ing inversely by the variance in the observed colours to 
obtain a mean best-fitting SED for each galax y. For the 
dust- c orrected colours we compare to the SEDs of lPoggianti 
( 1997) with 6 galaxy classes ranging from early-type galaxies 
(E/SO) to late-type spirals (Sc/Sd) supplemented with two 
starburst (SB) models, generated by superimposing a star- 
burst on a passively evolving system (see T98 and S2000 
for details). For our u ncorrected colours we use the SEDs of 
IColeman e t al. ( 1980), with 4 classes of SEDs from elliptical 
(E) to irregular (Im). 

In Fig. we present the results of the fitting process, 
showing the distribution of galaxy types (T) using each set 
of SEDs, with T = 1 representing the reddest SED. The dif- 
ferent classes are summarised in Table |5] Fig. Q shows the 
majority of the galaxies are types Sb or later, with very few 
early-type objects. There is a good cor respondence betw een 
the derived types (see Fig. |2Jl, with the lPoggiantil lll997l) de- 
rived types on average slightly bluer, possibly due to the in- 
clusion of dust corrections in the derivation of these types. 
The mean types in each case are T = 3.15 for the CWW 
SEDs (~ Sd) and T = 6.43 for the Poggianti SEDs (~ Sd). 
No galaxies appear redder than the model SEDs, though 
at the blue end, 3 of the galaxies are bluer than the bluest 
CWW SED, whilst 2 are bluer than the bluest SB2 Pog- 
gianti SED. This points to a large population of actively 
star-forming or starbursting galaxies in the sample, as ex- 
pected. 

One caveat is that these galaxy types (and hence fc- 
corrections) are based on the weighted average of (up to) 
four colours, and hence will not be sensitive to 'peculiar' 
colours in a single pair of filters where the other colours of 
that galaxy are 'normal'. This might occur (for example) in 
a galaxy experiencing a sudden burst of star-formation; the 
UV — B colours might appear very blue whilst the B — r' or 
r — i colours could be almost unchanged from their pre-burst 
values. We return to this issue in later sections. 
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Table 2. The correspondence between the different galaxy 
classes. 



SED set 



SO Sa Sb Sc Sd Im/SBl SB2 



Poggianti 1 
CWW 1 



6 



7 




Figure 1. The distribution of galaxy types (T) in our survey. 
UPPER: T = 1 corresponds to an E S EP, T = 8 corresponds 
to an SB2 SED. SEDs from IPoggiantl il997tV LOWER: T = 1 
corresponds to an E SEP, T = 4 corresponds to an Im SED. 
SEEs are from IColeman et al.l <1980l) . 
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3.3 Aperture Corrections 

In S2000, we saw that Ha derived SFRs were generally 
smaller than UV-based SFRs by factors of around 3-4. Po- 
tential (non-physical) explanations of this effect could be 
related to uncertainties in the flux calibration of the spec- 
tral data, or if the 2.7" diameter fibres used on the WYF- 
FOS spectrograph were sampling only a small fraction of the 
total galaxy light, and hence generating an aperture mis- 
match when compared to the FOCA magnitudes. In S2000, 
we found little significant trend in the ratio of Ha/UV light 
with redshift (and hence apparent size), as would be ex- 
pected if aperture corrections were not a significant source 
of uncertainty. However, the new B and r'-band data allow 
us to address any remaining uncertainties, as well as to ex- 
tend the number of objects in the survey with a calibrated 
Ha measure. 

For every FOCA B and r' counterpart, we measure 
from the imaging data the magnitude of each object in- 
side a radius of 2.7" (the WYFFOS fibre radius) using sex- 
tractor. We convolve the flux-calibrated spectrum with 
the filter+CCD response of the B/r' filters to generate cor- 
responding 'spectral magnitudes'. We find that the spec- 
tral B-r' colours typically agree with the optical 2.7" B-r' 
colours to within 10-20 per cent, indicating little wavelength- 
dependent uncertainty in the flux calibration. We then cal- 
culate aperture corrections by comparing the r' spectral 
magnitudes with the total r' magnitude (as measured by 
SEXTRACTOr), deriving correction factors by which we scale 
the Ha fluxes. We see minor evidence for a magnitude de- 
pendent effect, with some apparently brighter, and likely 
nearer with a greater apparent size, objects requiring larger 
correction factors. We exclude galaxies which require very 
large corrections from our sample and the future analyses 
in this paper. The median factor for the survey is 1.26, i.e. 
0.25 mag. For the small number of galaxies without any r' 
information this median correction is ap plied. This median 
factor compares well with that derived bv lTresse fc Maddoxl 
( 1998), who find a mean correction of 0.52 mag using V band 
magnitudes for the CFRS galaxies, but with smaller 1.75" 
slits, thus accounting for their slightly larger aperture cor- 
rections. 

The r'-band data also allow us to apply a 'flux cali- 
bration' to Ha fluxes which were taken during an observing 
run for which no standard stars were observed (detailed in 
T98), following a similar procedure as with the aperture cor- 
rections. This adds a further 17 objects to our Ha sample. 

These corrections assume that the Ha emission in the 
galaxy is a uniform distribution, and that the galaxy is not 
dominated by either nuclear or disc star-formation. Whilst 
these uncertainties may be significant for local nearby ob- 
jects (with a larger apparent size), by redshifts of 0.05 (i.e. 
for most of the galaxies in our sample) these effects are likely 
to be very small. 



Figure 2. Comparison of the galaxy types derived from the two 
SEP sets, with the CWW SEE types on the x-axis and the Pog- 
gianti SEE types on the y-axis, demonstrating the agreement be- 
tween the galaxy types derived from the two sets of SEDs, as well 
as before and after dust corrections have been applied. The stars 
show locations where the two types approximately correspond. 



3.4 UV/Optical colours 

One of the most puzzling features of the FOCA galaxies 
presented in our earlier papers was the presence of a sub- 
sample of galaxies with extremely blue UV — B colours, 
bluer than most typical starburst SEDs. The improved op- 
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Figure 3. Colour-redshift distributions for the unique optical 
counterpart cases only. UPPER: UV - B, LO WER: B - r The 
solid lines show the model predictions for the Pog giantl Jl997f) 
set of SEDs, the dotted lines are the colours derived from the 
CWW SED's (see text for details). The galaxies appear slightly 
bluer in UV — B compared to the model SEDs than for B — r; 
additionally in UV — B there is a small population of objects with 
colours bluer than the bluest model SED. The colours are in the 
AB system and are corrected for dust extinction. 



tical data that we have assembled now allows us to examine 
the SEDs of these objects across a wider wavelength range. 

With the new optical magnitudes, 8 and 5 per cent of 
the UV galaxies with a unique optical counterpart (slightly 
more when including multiple counterpart cases) are found 
to be bluer than the bluest Poggianti SED in UV — B and 
UV — g respectively (where the g data also includes the 
A1367 field). Before dust correction, these fractions are ~ 4 
and ~ 2 per cent respectively. At longer wavelengths, the 
galaxies have more 'normal' colours typical of the mean 
SED of the sample (see Fig- • As mentioned earlier, this 
might occur in galaxies experiencing a sudden burst of star- 
formation (see Section |2J . A second possibility is that that 
hot Wolf-Rayet (WR) s tars could be r espon sible for the UV 
excess as suggested by iBrown e t alJ teOOf J). We examined 
this suggestion in detail in JContini etld"N 2002h studying 
the co-added spectrum of the bluest objects in the sample, 
but could identify no WR spectral features in the spectrum 
of galaxies with extreme UV — B colours. 



4 STAR-FORMATION PROPERTIES OF THE 
SAMPLE 

We now turn to the main focus of this paper, continuing 
our investigations into the star-formation properties of the 
survey begun in S2000, where we examined the relationship 
between Ha and UV derived SFRs. We found that this rela- 
tionship possessed a scatter and non-unity best-fitting slope 
which could not be satisfactorily explained in scenarios as- 



Table 3. Default conversions between luminosity and SFR 

Tracer Luminosity for SFR = 1 M yr 

ergs - 1 A 1 ergs _1 Hz _1 
Ha 1.22 X 10 41 ergs" 1 — 

UV 2000A 5.75 x 10 39 7.78 x 10 27 

u' 3600 A 1.71 x 10 39 7.40 x 10 27 



suming constant or other simple SFHs. We now use our new 
data to investigate these effects with increased precision. 

4.1 Uncertainties in converting luminosities into 
star-formation rates 

Our first step is to construct a framework within which we 
can transform our measured luminosities for each of the star- 
formation diagnostics into SFRs. As has been covered ex- 
tensively in other studies, there are a remarkable number 
of free parameters in these models - the initial mass func- 
tion (IMF), the stellar metallicity, and the choice o f SED 
libraries, amongst many others (see lKennicuttl ll998. for a 
discussion of the issues). We derive conversion factors in a 
self-consi stent manner using the pegase- i i spec tral synthe- 
sis code dFioc fc Rocca-VolmerangeilT997lll99d) . The three 
star-formation diagnostics available in this study are 

• Ha luminosity (Lu a ), originating from re-processed ion- 
izing radiation at wavelengths A < 912 A produced by 
the most massive (> 10 Mq), short-lived (~ 20Myr), OB- 
type stars. Accordingly, Ha emission is a virtually instanta- 
neous star-formation measure, reaching a constant level after 
~ 10 Myr in a constant SFH. However, it is very sensitive 
to the form of the IMF due to the strong dependence on 
massive stars. Most calibrations assume case-B recombina- 
tion (a comprehensive trea tment of its use can be found in 
ICharlot fc LonghettillioOll) . 

• UV 2000 A continuum luminosity (Lxsv), originating 
from stars spanning a range of ages (and hence initial 
masses), including some post-main sequence contribution. 
Hence, any Ltjv to SFR calibration is dependent on the 
past history of star-formation, introducing a significant un- 
certainty when interpreting UV observations of star-forming 
galaxies. Dust corrections are also important in the UV, and 
a rel iable correction can be difficult to achieve (e.g. iBelll 
l2002ft . 

• u' 3600 A continuum luminosity (L u i), with a similar- 
physical origin as Ltjv- The principle disadvantage is that u' 
luminosities are contaminated by older, less massive stars to 
a greater extent than the UV-2000 A data (given the longer 
wavelength) , and this makes the calibration less certain than 
with shorter wavelength data. The advantage of using longer 
wavelengths is that dust corrections are smaller. 

Using pegase-ii, we investigate the dependence of the 
three diagnostics to the principle free parameters in the spec- 
tral synthesis codes, varying each to determine the spread 
in the conversion values. We investigate metallicities of 
Z = 0.0004, 0.0 04, 0.008 0.02 ( = Zq) and .05, together 
with t he IMFs of lSalpeterl ill955MScalol Jl99d) and lKroupal 
fcOOll) using mass ranges of 0.1 -> 100 M Q , 0.1 -> 120 M , 
0.1 -> 50 M Q , 1.0 -> 100 M and 5.0 -» 100 M . 

The uncertainties in the star-formation estimates due to 
these parameters are shown graphically in Fig. 2] as the dis- 
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Figure 4. The distribution of the conversion values for the Ho, 
UV and «' diagnostics (UPPER PANEL), and the distribution of 
the ratio of those values (LOWER PANEL), generated across the 
full range of metallicities, IMFs and mass ranges (see text for full 
details). The luminosities are those expected from star- formation 
rates of lMgyr -1 . The shaded histograms refer to conversions 
taken at 10 (heavy shading), 100 (light shading) and 1000 Myr 
(no shading) in a constant SFH. The figure demonstrates how 
the ratios involving Ha and UV (or it') evolve according to the 
time after the onset of star-formation, whereas the ratio of u'/UV 
has little time dependence. 

tribution of Lhq , Ltjv and L u / and the distribution of the ra- 
tios of Luv/Lhh, L u ' /Lhh and L u / /Ltjv for a constant SFR 
of lM yr _1 . There is a large spread in the conversion values 
for all diagnostics when considering the full range of IMFs 
and metallicities (around 1.5 orders of magnitude), with the 
importance of the timescale of recent star-formation for UV 
and u' diagnostics clear. However, the scatter when consider- 
ing the ratios of the diagnostics is much smaller, again with 
a strong time dependence. Though varying metallicities or 
IMFs can have a large effect on individual conversion values, 
the effects on the ratios is considerably smaller. Our default 
conversion values (listed in Table [3J assume solar metallic- 
ity, a Salpeter IMF with a mass range of 0.1 to 100 M©, and 
are taken 100 Myr into a constant SFH. 

4.2 Comparing star-formation diagnostics 

We present the relations between the different diagnostics 
in Fig. showing comparisons between dust-corrected Ha, 
UV and u' SFRs, overlaid with lines expected assuming con- 
stant star-formation and associated uncertainties in the po- 
sition of this line derived from the analysis of the last sec- 
tion. Also shown are the error-weighted least-squares best- 
fits to the data. The Hev-UV plot is an updated version of 
fig. 13 in S2000 with the new aperture and fc-corrections 
applied, together with some new measures and upper lim- 
its for Ha. Linear correlation (Pearson's r) and Spearman 
rank-order correlation statistics, together with the weighted 
least-square best-fit equations and associated x 2 i are pre- 



sented in Table 2] We do not plot the relations before dust 
correction, and instead list the same statistical data for these 
samples in Table 0] as a comparison and so the effect of our 
dust-correction prescription can be seen. 

These plots confirm several of the trends seen in S2000. 
Firstly, as in our previous studies we find correlations across 
three decades in luminosity between the different diagnos- 
tics. In all cases, the correlation coefficients are greater than 
0.90 (before and after dust correction). Secondly, the plots 
involving Ha-derived SFRs show a marked increase in scat- 
ter when compared to the it'-UV plot (see the x 2 figures for 
the best-fits in Table [IJ. The lines derived from the analysis 
of Section T4. 1 1 indicate that this scatter is unlikely to be pri- 
marily generated by varying IMFs and metallicities. Thirdly, 
the relation between Ha and UV-derived SFRs is complex, 
appearing luminosity (or SFR) dependent with the best- 
fit lines possessing non-unity slopes. At higher luminosities, 
the diagnostics agree well, whereas at lower luminosities the 
Ha diagnostic typically under-estimates the SFR when com- 
pared to both u' and UV continuum measure as shown by 
the slopes of the best-fitting equations in Table [I] The com- 
parison between the UV and u' diagnostics shows little lu- 
minosity effect with no offset from the default metallicity 
and IMF. Finally, Table 2] suggests that the scatter that we 
see is not a result of an inappropriate dust correction be- 
ing applied; the scatter in the relations decreases after our 
standard dust correction is performed. 

Further discrepancies from relations expected as a result 
of regular SFHs are also seen in panel D of Fig. Q>] which 
compares UV — B colours with the Ha equivalent widths. 
For a given Ha EW, many of the galaxies appear bluer than 
would be expected in simple SFHs. Again, simply varying 
the IMF does not appear able to reproduce the effect. 

As the UV and u'-band luminosities are measured at 
neighbouring wavelengths, they have a very similar depen- 
dence on the recent SFH (or the timescale over which the 
SFR varies) - see Fig. [I] The dominant cause of scatter in 
the UV — it' relation will be due to a combination of incor- 
rect dust-extinction corrections together with observational 
uncertainties (and the possibility of contamination by non 
star- forming galaxies), rather than varying star- formation 
timescales. The u'-UV relation for galaxies with detected 
Ha emission (i.e. known to be star-forming) is tighter than 
for galaxies selected without regard to Ha emission - indeed, 
many of the galaxies which show a large discrepancy between 
the UV and u' emission show no detectable Ha emission. 
Possibly these galaxies are not actively star-forming, with a 
substantial fraction of the UV and u' light arising from the 
presence of older and evolved, possibly post main-sequence, 
stellar populations. 

Though the luminosity (or SFR) dependent effect that 
we see in the plots involving Ha-emission could be due 
to a luminosity-dependent dust extinc tion relation (e.g. 
iHopkins et al 1 l200ll : ISullivan et afll200ll) , the correction wc 
applied in Section 13.21 already does a good job of remov- 
ing such effects in th e radio-UV relations for these galaxies 
JSullivan et alJl200ll) . making it unlikely the strength of the 
luminosity-dependent dust relation is under-estimated. 
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Figure 5. A comparison of A) UV2000A and Ha SFRs, B) u' 3600 A and Ha SFRs, C) UV2000A and u' 3600 A, and D) the UV-B 
colour and the Ha EW. All measures are dust corrected (and for Ha aperture corrected) as described in the text. In the UV2000 A — 
«' 3600 Aplot, galaxies possessing Ha emission are marked as filled circles, those without as empty circles. The SFRs are calculated at 
100 Myr into a constant SFH. The short-dashed lines denote an equality of SFRs and the long-dashed lines show typical uncertainties 
in the position of the short-dashed lines due to IMF/metallicity effects. Solid lines indicate weighted least-squares best-fits. In panel D), 
the overlaid lines correspond to the relation expected in a exponentially declining SFH for the three IMFs discussed in the text. 



Table 4. The results for the various statistical tests for Fig. [S] 
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4.3 The nature of star-formation in the sample 

The possibility of complex SF histories as been discussed in 
the past as an explanation for t he broad properties of late- 
type star- forming galaxies (e.g. ISearle et al.lll973l) . These 
studies have shown that periods of enhanced star-formation 
relative to a galaxy's history as a whole are able to reproduce 
well the blue colours of local dwarf galaxies, as well as the 
scatter that is observed in the colours of these galaxies. In 
S2000, we introduced the possibility of complex SF histories 
as a possible explanation for the discrepancies between Ha 
and UV luminosities. The results of the last section based 
on additional u-band information, an extended optical wave- 
length coverage, improved fc-corrections, and aperture cor- 
rections on the Ha fluxes, confirm these findings. In S2000, 
we introduced the concept of varying SFHs as a possible ex- 
planation for our dataset, and we now investigate this fur- 
ther. 

The key is the time-scale upon which different diag- 
nostics of star-formation trace changes in a galaxy's SFH. 
The Ha luminosity depends only on the most massive and 
short-lived stars, and the point at which new stars are born 
at the same rate as older stars die occurs after only a few 
million years. By contrast, as UV/ti' continuum measures 
have a significant contribution from older and longer-lived 
stars, it takes 100 to 1000 Myr to reach the stage at which 
the birth-rate of the stars which generate UV/w' emission is 
the same as their death-rate. A burst or increased period of 
star-formation superimposed on an otherwise regular SFH 
therefore affects the different diagnostic plots in different 
ways. 

For the Ha versus UV (or u') plot, a star-formation 
event will move a galaxy rapidly up and right on the diag- 
nostic plot with a short period in which the observed Ha 
SFR is higher than the observed UV SFR as the UV light 
catches up with the Ha light. As the burst subsequently dies 
away, the Ha rapidly decreases, whilst the UV luminosity is 
retained due to the contribution from older stars. The galaxy 
describes a loop in Ha-UV (or Ha-it') space (see panel A 
in Fig. ||)J . The size of the loop depends on the parameters 
of the burst: stronger (more massive) and/or shorter bursts 
produce larger loops. By contrast, the effect on the UV-ti' 
relation is small as these diagnostics have a similar depen- 
dence on the SFH. For the UV — B versus Ha EW relation, 
the burst rapidly increases the Ha EW and generates a bluer 
UV — B colour. As the burst dies away, the Ha returns to 
pre-burst levels, whilst the UV — B colour remains blue (see 
panel B in Fig. |HJ . 

This explanation is supported by our observations and 
also predicts trends we should see in our spectral diagnos- 
tics indicative of the age of the stellar population. In galax- 
ies which are either at the peak of a starburst or which are 
undergoing a regular star-forming process (those galaxies in 
which SFRs derived from Lh and Luv agree well) , we would 
expect small values of the Balmer decrement and high values 
for the Ha EW due to the presence of a young stellar popu- 
lation. For galaxies which are in the later stages of a burst of 
star- formation, where we have excess in UV/it' SFRs com- 
pared to Ha, we would expect larger Balmer breaks and 
weaker Ha EWs as the burst population ages. 

We examine these parameters in Fig. EI showing the dis- 
tribution of the Balmer break (D4000) and Ha EW values as 



0.6 
0.4 
0.2 



E I 



10 50 100 

Ha equivalent width (A) 



0.6 
0.4 
0.2 


-0.2 



• . - *" . • 



1.2 1.4 1.6 1.1 

D4000 



Figure 7. The distribution of the residuals from the constant 
SFR line on the HagpR— UVgpR plane, as a function of the Ha 
EW (upper) and D4000 (lower) for the FOCA galaxies, shown 
as solid circles. The horizontal dashed lines show the position of 
constant star-forming galaxies. Overlaid stars show the evolution 
of these parameters during example starburst events as in Fig. 
Again, the starburst events are able to reproduce the scatter and 
offsets seen. 



a function of the position of a galaxy on the HasFR - UVsfr 
diagram. As expected, we see some evidence that those sys- 
tems with the largest residuals from constant SFR scenarios 
demonstrate features indicative of older stellar populations 
(larger values of D4000 and smaller Ha EWs). 

Irregular SFHs clearly help to explain the discrepan- 
cies that we see in the relations exploring the different star- 
formation diagnostics, in particular the 'excess' SFRs de- 
rived from UV and u' measures when compared to Ha, and 
the luminosity dependence of this result. We investigate this 
hypothesis more quantitatively in Section |S] 



5 ANALYSIS 

We have reviewed the qualitative evidence that many of the 
low to intermediate luminosity galaxies in this data set do 
not possess simple constant or smoothly declining SFHs. 
This is an important finding in the context of the interpre- 
tation of SFR measures from large galaxy redshift surveys, 
our primary motivation in this section is to quantitatively 
confirm our findings. With sample photometry that ranges 
from 2000 A to ~ 7000 A (plus spectral coverage at optical 
wavelengths), we are able to investigate any SFHs that are 
more consistent with the combined dataset, via a simple x 2 
modelling technique. We do this by examining the form of 
the recent SFH in our galaxies, and then discuss the impli- 
cations for utilising the different star-formation diagnostics 
for the galaxy population as a whole. 
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Figure 6. The effect of an enhanced star-formation event on A) the Lft a -Ltjv relation, and B) UV — B versus Ho EW. Two example 
bursts are shown: one involves 10 per cent of the galaxy mass in an exponentially decaying burst of e-folding time 45 Myr, the second 
30 per cent mass with a time of 20 Myr, generated using the PEGASE-II code, a standard IMF, and solar metallicity. The enhanced 
star-formation events are able to reproduce the scatter and offsets from constant star-formation scenarios seen in the data. 



5.1 Evidence for complex star- formation 
histories? 

For each galaxy we have available the following (dust- 
corrected) measures: i) The Ha luminosity, ii) The UV con- 
tinuum luminosities at 2000 and 3600 A, iii) The strength of 
the Balmer Break, D4000, iv) The Ha EW, v) The UV - B 
and B — r' colours. 81 of our galaxies fit these criteria. The 
amount of extinction in the models is not left as a free pa- 
rameter, but is constrained observationally using the Ha/H/3 
Balmer decrement as described earlier. Our goal is to ex- 
amine the form of the most recent star-formation event for 
which this is an excellent dataset (however it is not able 
to provide constraints on the SFH in the distant past). In 
what follows, we assume a Salpeter IMF with mass limits at 
0.1 Mq and 100 Mq, solar metallicity and again utilise the 
PEGASE-n population synthesis code. 

We begin by considering two simple SFHs, one repre- 
senting a constant SFH and the other an exponentially de- 
clining SFH according to 



SFR(i) = M (M yr 



SFR(t) = ^exp(-l) (Moyr- 1 ) 

T T 



(3) 
(4) 



In equation Q, Mo is the SFR and the stellar mass formed 
at time t is Mo x t; in equation Q, r is the time constant 
of the SFH (in Myr) and Mo is the total mass formed at 
time t — oo (i.e j Mo/t exp(— t/r)dt = Mo). Small values of 
t (e.g. t < 5 Myr) approximate instantaneous bursts, large 
values approximate continuous SFHs. 

We construct a grid of galaxy spectra corresponding to 
the different SFHs as a function of time, with a range of r 
from 500 to 5000 Myr for the exponential model and models 
ages of up to 12 Gyr. For each galaxy, we shift the model 



spectra to the observed redshift, and calculate predictions 
for each of the observed parameters listed above including 
nebular emission according to the prescription in pegase-ii. 
We normalise each model (i.e. adjust Mo) so that the pre- 
dicted and observed Ha emission are identical. We then cal- 
culate a v 2 statistic for each model galaxy by comparing the 
observed and predicted parameters, and find the combina- 
tion of parameters which minimise this y 2 - The probability 
of each y 2 is then calculated using the incomplete gamma 
function for the appropriate number of degrees of freedom 
in each fit. 

The key result is that the 'best-fits' for 46/81 (57 per 
cent) of the galaxies in the constant SFH scenarios and 37/81 
(46 per cent) of galaxies in the exponential SFH scenarios are 
rejected at greater than the nominal 99.9 per cent confidence 
level - i.e for just under half the sample, the model SEDs 
provide poor representations of the true physical scenario 
over the full range of t and r. The principle indicators of 
the poor fits are the u'/UV luminosities and the HaEWs, 
sensitive to the ratio of ongoing to past star-formation. 

Consequently, for the 46 galaxies which are poorly fit 
by simple SFHs, we extend these models by relaxing the as- 
sumption of 'smooth' SFHs and include the possibility of a 
burst of star-formation superimposed on top of the under- 
lying galaxy SFH, represented by 



SFR(t) 
SFR(t) 



(1 ~ h)Mp t 
exp( ) if t < t b 

(1 - f b )M , t 
•i '- exp( ) 

To T 



•a 
f b M 

Tb 



exp(- 



tb — t , 

Tb ' 



if t >= t b 



(•») 



with the starburst commencing at time tb and the total 
mass formed Mo, a fraction formed during the starburst. 
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Table 5. The median best-fitting parameters for the form of the 
bursts of star-formation 

Parameter Median value 

7 b 23T 

t b 125 

fb 0.09 

LOG M stcllar 9.91 



We ignore times t < t b (simulated above), and assume 
tb = 8000 Myr and r g — 3.0 Gyr. The results are not sen- 
sitive to the choice of these two parameters which define the 
underlying SED. 

We fit for the combination of parameters which min- 
imise the x 2 as before, with stellar masses estimated by in- 
tegrating the SFH to the best-fitting t. The results of the 
fits using this framework is quite different - all but 6/46 of 
the galaxies which were poorly fit by regular SFHs now have 
acceptable fits. Of these, only 3 have fits which are excluded 
with a very high level of confidence with the 3 others 'bor- 
derline' fits. The median best-fitting parameters are given 
in Table |^| Though these parameters do vary widely when 
considering any one object, they suggest that the most com- 
mon starbursts are relatively short duration events, involv- 
ing around 8-10 per cent of the galaxy mass, and that we 
view these galaxies not at the peak of a particular starburst 
(where Ha and UV derived SFRs will be approximately 
equal), but instead some way into a burst's lifetime; it is 
at these times that the Ha and UV or u' luminosities are 
the most discrepant. 

Our conclusion from this modelling work is that whilst 
regular SFHs can provide an adequate explanation for ap- 
proximately half of the sample, they provide poor represen- 
tations of the remaining galaxies. By relaxing the assump- 
tion of a simple SFH, and instead allowing the recent star- 
formation in the galaxies to evolve according to a simple 
burst structure, these inconsistencies appear to be resolved. 

5.2 Comparisons with other samples 

The finding that the SFHs of a substantial fraction of the 
galaxies in our dataset appear irregular has important im- 
plications for surveys that measure SFRs in galaxies via 
UV measures, which we discuss in the next section. Prior 
to that analysis, we first examine the results of this sur- 
vey with other redshift surveys as a consistency check to 
confirm that the UV-selected galaxy properties make sense 
within the broader galaxy population. 

Comparison of the results of this survey with other 
samples of UV-selected galaxies at low-redshift are not yet 
possible due to the lack of low-redshift UV observations (a 
situation soon to be rectified via the GALEX experiment), 
and [/-band observations must suffice. However, at higher 
redshift a comparison is easier as (7-band fluxes are shifted 
into optical bandpasses. We compare here with five different 
samples. The first, at z ~ 0.2 5 and z ~ 0.4, is the em ission- 
line (Ha) selected sample of [Hippclcin_£t_alJ (|20Qj) . Next 
are two high-redshift samnlesrf lGdazebrooket^Ll (|l999T ) 
(z ~ 1), who observe Ha fl uxes for a samp le of /-band 
selected CFRS galaxies, and lErb et al.l J2003lh who select 
z ~ 2.2 galaxies using a photometric col our technique. Fi- 
nally, we show the two z ~ samples of iBell fc Kennicuttl 
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Figure 8. A comparison of the FOCA < z < 0.4 galaxies 
(data same as Fig. |SJ with those from other sa mples at varying 
redshi fts: Open circl es: z ~ 0.25 and z ~ .4; iHippelein et al] 
i2003ft. Stars: z ~ 1: iGlazebrook et al.Ul99SD. Squares: z ~ 2.2 ; 
lErb et alJ<2003l) . Square crosses: z ~ 0:lBell fc Kennicuttl <200ll) . 
Triangular crosses: z ~ 0: lBuat et al"l<2002|y The solid line shows 
the best- fit to the FOCA sample. All data are dust-corrected using 
the lCalzetti et all J200fj) prescription. 



<200ll) and iBuat et alJ <2002l) . For all these samples, we 
take the published Ha and UV- c ontinu um fluxes, correct 
for dust using the ICalzetti et al] 1120001) prescription if re- 
quired, and convert to SFRs using our cosmological model 
and the PEGASE-n conversions as appropriate. The samples 
are then plotted with the FOCA data in Fig. |H] 

In high- luminosity systems (SFR > 5MQyr _1 ), the 
only systems probed by the high-redshift studies, the differ- 
ent samples agree well, with some evidence for higher Ha- 
derived SFRs in these systems. Whilst this could be caused 
by under-estimated dust corrections, our alternative hypoth- 
esis of irregular SFHs can also explain the observations: high 
Ha luminosity systems are likely near a peak in SFR, where 
models predict that Ha and UV luminosities are at least 
equal, or should even show an Ha excess (as Ha light in- 
creases in a new starburst more rapidly than UV light). 

Below SFRs of ~ 5Mq yr _1 , the properties of the dif- 
ferent samples begin to diverge. As noted in Section 14.21 
the UV-selected sample shows a general excess in the UV- 
derived SFRs. However, the two z ~ samples demonstrate 
SFRs that agree well, while the Ha-selected sample shows 
a larger scatter in the derived SFRs, with an Ha excess at 
the faint end. Once again, a systematic under-estimation in 
the dust extinction corrections can explain the Ha-excess 
at the faint end, but this does not explain the UV-excess 
galaxies. Again, all of the results can be consistently inter- 
preted within a framework of non-regular SFHs coupled with 
the various survey selection criteria: Ha-selected galaxies are 
likely to be located at a phase in their SFH where they are 
either near a peak of a starburst (in galaxies with varying 
SFHs) and the Ha-derived SFR is greater than that from the 
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UV, or the galaxies will have regular SFHs and the Ha and 
UV SFRs agree well. The z ~ galaxies, with a general op- 
tical selection criteria, are likely to represent star-formation 
across normal Hubble types and are less likely to possess 
the irregular SFHs found in the FOCA sample which are 
selected by their UV light. We examine these ideas further 
in Section f5.3l 

5.3 Implications for measuring cosmic 

star-formation for flux-limited redshift 
surveys 

The finding that a fraction of the galaxies in our sample 
do not possess simple SFHs has implications for flux-limited 
redshift surveys such as this. A galaxy that evolves with 
an intermittent SFH will clearly brighten and dim over the 
course of its history. We demonstrate this in Fig. [5] where 
we show the UV evolution of a typical galaxy in our sur- 
vey and the input SFR(t) required to produce it, together 
with the flux limit of the FOCA survey (m uv = 20.75 in 
the AB system). We also plot the flux limit and survey pa- 
rameters of the z ~ 0. 25 and z ~ 0.4 Ha-selected survey of 
IHippelein et~ai] J2003t) . 

The bursts of star-formation brighten the model galaxy 
above the UV flux limit into the detectable magnitude range. 
Without this boost, a galaxy would otherwise not be de- 
tected in the UV unless located at a low redshift. This bias 
in itself is not serious - it is one of the goals of this red- 
shift survey to measure an integrated star-formation den- 
sity, including that fraction of star-formation that occurs in 
a 'burst mode' rather than a 'continuous mode'. However, 
a more important bias arises due to the time required for 
the UV light to die away after a burst has completed, thus 
leading to higher measured SFRs in the UV or u' than might 
be obtained via alternative diagnostic measures, and higher 
measured SFRs than the 'true' SFR in the galaxy. 

We illustrate this using a model taking into account the 
flux limit of the FOCA survey and Ha-selected survey. We 
model a galaxy's SFH by superimposing a series of bursts 
onto an exponentially declining SFH between the ages in 
the galaxy's history that correspond to * = 0.5 and z — 
in our cosmological model. The variations UV(t) and Ha(£) 
are calculated, and using the mapping of time — * z, the ap- 
parent magnitude m uv (t) can then be computed from the 
luminosity distance and an appropriate fc-correction calcu- 
lated from the synthetic spectrum. From the flux limit of 
the two surveys, the visibility of the galaxy in each survey 
can be found. 

Fig.|5]shows this model. There are points over the model 
galaxy's evolution when, although the instantaneous SFR 
of the galaxy lies below the flux limit of the UV survey - 
i.e. if this SFR were converted to a UV luminosity using a 
simple constant-SFH assumption the galaxy would not be 
detected - the galaxy remains in the UV-selected survey 
due to the slower decline of the UV light. As expected, at 
higher redshift a galaxy is preferentially picked out if it is 
undergoing a burst of star-formation and additionally if the 
galaxy is near the peak of the particular starburst. At lower 
redshifts, we are able to view the galaxy at later times into 
a particular starburst event, and only at the lowest redshifts 
does the underlying (smoothly declining) galaxy SED come 
into the FOCA survey. 
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Figu re 9. The vi sibility in the FOCA UV-selected and 
IHippelein et alj i2003t) emission-line selected surveys of a model 
galaxy undergoing a series of starbursts between z = 0.5 and 
z = 0. The continuous line denotes the UV luminosity; the dashed 
line the input SFR, (the axes are scaled so that the UV luminosity 
and input SFR agree according to the conversion values in Ta- 
ble 0- The filled circles show the periods during which a galaxy 
falls within the selection criteria of a FOCA-like survey, and the 
two overlaid dotted curves show the magnitude limit of the FOCA 
survey (upper curve) a nd the sensitivity limit of the Ho survey of 
IHippelein et all 120031) (lower curve) . The vertica l shaded areas at 
~ 0.25 and z ~ 0.4 show the redshift sensitivity of lHippelein et all 

The effect of these forms of SFHs on the Ha-UV plane 
in a UV-selected survey can be seen in Fig. If Ul We run 200 
simulations, generating SFHs containing a random number 
of bursts (from zero to three) of varying mass and duration 
superimposed on underlying histories with r from 0.750 to 
6Gyr between z = — 0.5. The Ha and UV luminosities 
at 1 Myr intervals are calculated and given a random error 
distribution identical to that in the survey. The galaxy is 
recorded (number-weighted by the co-moving volume) if it 
meets the selection criteria of the UV survey, and then exam- 
ple survey datasets drawn randomly from these results. The 
distribution of 18000 simulated galaxies in Ha-UV space 
can then be generated. An excellent match between the ob- 
served galaxy distribution and the simulated galaxies can be 
seen. 

Estimating how 'biased' a determination of the SFR 
derived from a UV continuum measurement in a partic- 
ular galaxy is (i.e. how much this measurement is over- 
estimated due to recent star-formation activity) clearly re- 
quires a knowledge of the precise form and age of the par- 
ticular star-formation events. Though these parameters can, 
in principle, be estimated using techniques such as those in 
Section f5. II this can only be done reliably for a sub-sample 
of the galaxies - and even then there is a considerable uncer- 
tainty in the derived parameters due to the lack of infra-red 
photometry (particularly the total stellar mass). 
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tions given above are only applicable to individual galaxies, 
not to the integrated luminosity density o f the survey as a 
whole, as the standard luminosity function ISchechterl jl976fl 
parameters (a, M,/L» and 0*) will be affected in different 
ways by correcting for this effect. This correction tends to 
make the SFR in fainter galaxies lower, and hence the ap- 
propriate luminosity function slope will be flatter, or a will 
become less negative, likely leading to a slight decrease in 
the calculated integrated star-formation density. The lumi- 
nosity functions and light densities of this sample will be 
addressed in a forthcoming paper (Treyer, Sullivan & Ellis, 
in preparation). 

The conclusion from this study is that integrated lu- 
minosity densities are a poor guide to the complex physical 
processes at play in these systems. Star-formation densities 
derived from these measures need to be carefully calibrated, 
and the effects of different biases at work more fully under- 
stood, before the results of flux-limited surveys can be fairly 
compared. 



Figure 10. An example simulation of the galaxies in our survey 
generated by superimposing starbursts of various masses and du- 
rations on a selection of underlying galaxy SFHs. By sampling a 
galaxy's history at random redshifts (weighted by the co-moving 
volume) during which the galaxy's apparent UV and Ha lumi- 
nosities allow detection in the FOCA survey, distributions in the 
Ljj q - Ljjv plane (panel A from Fig. |5J can be generated and 
compared to the observations. Simulated galaxies are denoted by 
small dots, observed galaxies by open circles. The solid line de- 
notes equality of SFRs from Table 151 



Instead, we can estimate the effect of the recent SFH in 
a purely statistical manner based on the entire sample. We 
can use the plots of Fig.|S]and the statistical tests performed 
on them to derive relations between UV or u' luminosities 
and the instantaneous (Ha-derived) SFR. By fitting the data 
across the range of SFRs probed in this study, we find: 



log(L 



uvj 



SFR 



= 40.03 + 0.75 log(SFR Ho 



(- 



JUV 



1.07 x 10 40 
for the UV luminosities, and 



) 1 ' 33 Moyr" 1 



(6) 



log(W) 
-> SFR 



39.53 + 0.74 log (SFRhci) 

/ Luv xl.35 



'3.39 x 10 39 ' 



M 



©yr 



(7) 



for the it' luminosities, where the luminosities are dust- 
corrected and measured in ergs -1 A . This correction is 
valid over SFRs of approximately 0.1 to 100 Moyr -1 for 
galaxies in this survey selected at 2000 A. Whilst similar 
relational forms likely hold for other redshift surveys, the 
exact form will depend on the selection criteria and hence 
the number of galaxies with 'normal' SFHs admitted into 
the survey. Even for other UV-selected surveys, selection at 
a shorter wavelength will be less affected by the effects of 
starbursts as the UV light dies away more quickly, whilst 
longer wavelengths could suffer considerable bias. 

Understanding the precise impact of this bias on analy- 
ses such as the Madau plot is a complex problem. The correc- 



6 CONCLUSIONS 

In this paper, we have presented the results of panoramic 
wide-field optical imaging of a UV-selected galaxy redshift 
survey, with the aim of further investigating the nature of 
star-formation in local z < 0.4 star-forming and starburst 
galaxies. We have found the following: 

(i) New u Br' photometry, supplemented by DPOSS 
data, have allowed us to derive aperture corrections for our 
Ha fluxes, as well as more reliable SEDs for all the sample 
galaxies. There remains a small fraction 8 per cent) of 
galaxies exhibiting UV — B colours bluer than most star- 
burst models, whilst possessing colours typical of standard 
SEDs at longer wavelength. 

(ii) We have investigated the dependence of the Ha, UV 
and it' luminosity to SFR conversion factors on the IMF and 
mass ranges, the stellar metallicity and the time since onset 
of star formation, assuming constant SFHs. Though varying 
these parameters can have a large effect on individual con- 
version values, the effect is smaller when considering ratios 
of these values (e.g. UV/Ha). 

(iii) T aking advantag e of ou r new dataset, we update the 
study of ISullivan et all (|200CT) and compare SFRs derived 
from UV or u' and Ha luminosities. Assuming simple SFHs, 
we show the scatter and non-unity best-fitting slopes ob- 
served are unlikely to be primarily generated by varying the 
above parameters. 

(iv) We show that models including a burst or increased 
period of star-formation superimposed on an otherwise 
smooth underlying SFH provide much better fits to the data 
set of ~50% of the galaxies (Ha, UV and u' luminosities, 
Balmer Break, Ha EW and colours) than a smooth SFH 
alone. 

(v) Such burst modes of star formation lead to an overes- 
timate of SFRs derived from UV luminosities in lower lumi- 
nosity systems, as UV light from less massive stars will still 
be present after a burst has died away. We propose a simple 
statistical correction for UV-selected surveys based on the 
'true' SFRs as derived from Ha luminosities which originate 
only from the most massive stars. 
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